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Perspectives : 


Bibliographie des techniques de programmation : 


Les algorithmes génétiques sont des algorithmes de recherche inspirés des 
mécanismes de l’évolution naturelle des êtres vivants et de la génétique. 


Ces algorithmes sont basés sur des opérateurs qui simulent l’évolution d’une 
population. Ils ont été utilisés avec succès dans divers domaines tels que 
l’optimisation d’emploi du temps, de contrôle de systèmes industriels, de réseaux 
d’antennes ou encore de commutateurs optiques adiabatiques. Ils ont aussi été 
utilisés dans le domaine de la mécanique, notamment pour l’optimisation de 
dimensionnement de structures composites stratifiées. 


Dans ce contexte, il s’avere important de s’interroger sur l’utilité des 
algorithmes génétiques ainsi que la méthode de monte Carlo dans |’ optimisation 
des structures dans le domaine du Génie civil. 


Pour étayer profondément ce rapport, nous allons voir successivement le 
principe général d’un algorithme génétique, son organigramme, son role et ses 
avantages enchainé par son application comparée a la méthode de Monte Carlo 
aux structures tridimensionnelles en acier, avant de mettre l’accent sur une 
application nette des algorithmes génétiques a une poutre-console en béton armé 
a inertie constante. 


Un algorithme génétique recherche le ou les extrema d’une fonction définie sur un espace de 
données. Pour l’utiliser, on doit disposer des cing éléments suivants : 

v Un principe de codage de l’élément de population. Cette étape associe à chacun des points 
de l’espace d’état une structure de données. Elle se place généralement après une phase de 
modélisation mathématique du problème traité. Le choix du codage des données 
conditionne le succès des algorithmes génétiques. Les codages binaires ont été très 
employés à l’origine. Les codages réels sont désormais largement utilisés, notamment dans 
les domaines applicatifs, pour l’optimisation de problèmes à variables continues. 


v Un mécanisme de génération de la population initiale. Ce mécanisme doit être capable de 
produire une population d'individus non homogène qui servira de base pour les générations 
futures. Le choix de la population initiale est important car 1l peut rendre plus ou moins 
rapide la convergence vers l’optimum global. Dans le cas où l’on ne connaît rien du 


probleme a résoudre, il est essentiel que la population initiale soit répartie sur tout le 
domaine de recherche. 


v Une fonction à optimiser. Celle-ci prend ses valeurs dans R+ et est appelée fitness ou 
fonction d’évaluation de l’individu. Celle-ci est utilisée pour sélectionner et reproduire les 
meilleurs individus de la population. 


Des opérateurs permettant de diversifier la population au cours des générations et 


d’explorer l’espace d’état. L'opérateur de croisement recompose les gènes d’individus existant dans 
la population, l’opérateur de mutation a pour but de garantir l’exploration de l’espace d’état. 


v Des parametres de dimensionnement : taille de la population, nombre total de 
générations ou critère d’arrêt, probabilités d’ application des opérateurs de croisement et de mutation. 
Le principe général du fonctionnement d’un algorithme génétique est représenté sur la 
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Figure 1 : Principe général des algorithmes génétiques. 


On commence par engendrer une population d’individus de façon aléatoire. 
Pour passer d’une génération k a la génération k +1, les trois opérations suivantes 
sont répétées pour tous les éléments de la population k. Des couples de parents 
PI et P2 sont sélectionnés en fonction de leurs adaptations. L’opérateur de 
croisement leur est appliqué avec une probabilité Pc (généralement autour de 0.6) 
et engendre des couples d’enfants Cl et C2. D’autres éléments P sont 
sélectionnés en fonctionde leur adaptation. L’operateur de mutation leur est 
appliqué avec la probabilité Pm (Pm est généralement très inférieur à Pc) et 
engendre des individus mutés P. Les enfants (C1,C2) et les individus mutés P 
sont ensuite évalués avant insertion dans la nouvelle population (la Figure 1 ) 
présente le cas ou les enfants et les individusmutés remplacent les parents. 
Différents critères d’arrét de l’algorithme peuvent être choisis : 
— Le nombre de générations que l’on souhaite exécuter peut-être fixé a 


priori. C’est ce que l’on est tenté de faire lorsque l’on doit trouver une 
solution dans un temps limité. 


— L’algorithme peut être arrêté lorsque la population n’évolue plus ou 
plus suffisamment rapidement. 
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Figure 2 : Organigramme de la méthode des AG 





Les AG ont trois (03) avantages : 
— l’optimisation de structure composite comporte de nombreux optima locaux. 
— les AG ne donnent pas une solution unique. En effet, une famille de 


solutions est obtenue, ce qui permet d’atteindre plusieurs optima. 
Cette particularité a pour conséquence de conduire à une optimisation riche 
en solutions potentielles. Elle permet également d’intégrer des contraintes en 
post- traitement qui ne peuvent pas être prises en compte directement dans 
PAG. 
— les AG permettent aussi bien l’optimisation de problèmes dépendant 
de variables discrètes que de variables continues. 






Pour montrer l’intérêt de la méthodologie proposée, nous allons adapter la 


stratégie qui consiste à évaluer la même structure par deux méthodes différentes : 
1. la première suivant l’approche d'optimisation de Monte-Carlo, 
2. la deuxième avec l’outil qui correspond à une optimisation par la 


méthode des algorithmes génétiques. 


Cette façon de procéder nous permettra de surmonter les deux grandes difficultés 
rencontrées le plus souvent par les programmeurs dans la résolution de problèmes 
physiques ou réels en l'occurrence la validation des résultats numériques obtenus 
et la performance de la méthode de résolution choisie . Nous considérons la 
structure métallique suivante (Figure 3): 
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Figure 3 : Structure en 3D avec des sections circulaires 


- Les connexions entre les éléments (n = 4) sont parfaitement intégrées (rigides) 


- Nous adoptons comme variables d'optimisation, les sections S1 et S2 qui 
désignent respectivement les sections des poteaux et celles des poutres ainsi 
que L1,L2 et H qui représentent respectivement la longueur, la largeur et la 
hauteur des semelles de fondation. 

- La limite élastique de l’acier est fy=235 MPa, 


- Le tassement différentiel admissible est égal a 0.1mm 
- La profondeur de fondation est fixée a 2m. 


- L'enveloppe financière réservée à ce projet est de: 124,5 à 130,2 DH. 


- Le chargement vertical est réparti sur les 4 poutres, tandis que le 
chargement horizontal est distribué partiellement sur un des 4 poteaux, 


Pour l’estimation du projet, nous avons adopté des barèmes de prix moyens 
quisont applicables dans le marché du bâtiment, soit : 
Tableau 4 : Barèmes de prix moyens de matériaux 


Désignation 


Barèmes (DH/tonne) 


Acier de 
Structure Acier 
de ferraillage 


Béton Armé 





La méthode de Monte-Carlo (en fait, toute une gamme de méthodes [Fishman, 1995] 
consiste à rechercher l’optimum d’une fonction, par une heuristique d'essais erreur, en 
générant une suite aléatoire de solutions suivant une loi statistique (en général uniforme). 
C'est la technique la plus simple à mettre en œuvre lorsque la taille de l'espace de recherche 
dépasse les capacités d'énumération disponibles. Cette méthode présente de bonnes 
capacités d'exploration mais aucune capacité d'exploitation. Elle nécessite un nombre 
d'évaluations très important [Yang,2006]. Pour pallier ce manque d'efficience, on 
privilégie les méta-heuristiques pour guider la recherche et acquérir de réelles capacités 
d'exploitation. Deux classes de méta- heuristiques sont présentées de façon succincte dans la 
suite : les méthodes de recherche par voisinages et les algorithmes évolutionnaires. 





- Écart local: représente la précision de Monte-Carlo. 

- Ecart minimum: est le meilleur écart de la norme euclidienne à l'itérationactuellei 
de MC. 

- Ecart actuel: est l'écart min de l'itération précédente. Si c'est la première 
itération, sa valeur est prédéfinie dans le code source et est égale à 10* (c'est 


le plus grand écart connu actuellement pour l'algorithme de Monte Carlo). 
- nit: représente le nombre d'itérations de Monte-Carlo. 
- Resumpt: représente le nombre de reprise de Monte Carlo. 
Les résultats obtenus par la méthode de Monte Carlo sont présentés dans leTableau 


1: 


Tableau 1 : L'optimisation Monte-Carlo 


Abaqus 
Numet Calcula |4 4 4 
-tion 


Par_nomgene | X57320 S1=35 S1=461.00 | S1=459.48 
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(mm2)and u S2=30 Length=758. | Length 


(mm) i 0 2Width -AIA 


Support —1000 
> Width =500 
Parameters 
variation ( % 
) 


Response 
values 
(MPa) 


Optimise S 1=461.00 S 1=459.48 S 1=464.64 
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parameter Length=75 8. Length=747. Length=742. 

s (mm2) 2 4 5- 
Width=363.3 | Width Width =352.4 
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4.22e™ cost 
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O1 

O2 
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27250.5diff 
“cost=846.21 


2 
4 
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MC 
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Le problème de la conception globale des structures en acier peut, à notre 





avis,étre globalement posé comme un problème d'optimisation pour minimiser la 


fonction objectif ou le critère du coût global (CG) de la structure conformément à 
la condition ou contraintes de l'Eurocode. 


Le (CG) est une fonction d'un certain nombre de variables telles que les trois 
variables: I, Xa et Xn et seront expliquées explicitement dans les sections 
suivantes. 


Ainsi, nous pouvons formuler implicitement le problème d'optimisation 


comme suit: Min CG (I, Xa, Xn) 


ou: CG: le coût global de production de la structure. 
I : vecteur des caractéristiques dimensionnelles des 
barres.Xa : vecteur nature des appuis. 
Xn: vecteur nature des nœuds. 
Le coût global (CG) d’une construction métallique (superstructure et 
fondations )peut s'écrire: 
(CG) = (CS) + (CF) 
(CS) : Coût Superstructure en acier. 
(CF) : Coût Fondations en béton 
armé. 
Le coût (CS) est lui-même composé des coûts suivants : 
(CS) = (Mat) + (Fab) + (Mon) 
(Mat) : Coût des matériaux des profils et 
assemblages,(Fab : Coût de la 
fabrication en atelier. 


(Mon): Coût de montage des différents éléments sur site. 
Le coût (CF) est composé des deux entités élémentaires 


suivantes : (CF) = (Ter) + (PrF) 


(Ter) : Coût des terrassements. 
(PrF) : Coût de production des fondations. Il intègre les coûts des matériaux 
etles coûts d'exécution. 





- Concernant l'approche de Monte Carlo, 1l faut noter que chaque itération de Monte Carlo 
nécessite 34 = 81 tirages (donc 243 calculs pour cette structure). 

L'exploitation des résultats de l'optimisation de Monte Carlo montre que la solution 
optimale est celle qui correspond au minimum de min-gap (ici égal à 2428.8). Par 
conséquent, la solution est le vecteur (S1 = 464.64, S2 = 181.96, longueur = 742.5, largeur 
= 352.4), 

L'écart local n'étant pas atteint, il s'agit donc d'un optimum local. Cependant, cette solution 
correspond à des contraintes d’acier raisonnables et au coût minimal (ici égal à 534,7 
DH). 

Concernant l'approche des algorithmes génétiques, 1l est notoirement connu que les 
algorithmes génétiques utilisent souvent des vecteurs et des matrices pour modéliser les 
problèmes d'ingénierie. 

Par conséquent, notre exemple d'application a été implémenté avec Java et utilisant un 
code de calcul basé sur des éléments finis ABAQUS. 

En suivant scrupuleusement les étapes du paragraphe ci-dessus et en appliquant les valeurs 
0,6 et 0,006 respectivement pour les probabilités de croisement et de mutation, le meilleur 
vecteur de solution trouvé par la méthode des algorithmes génétiques après 27 itérations 
est: 


V=(21, 10, 31, 18) ce qui correspond à 40DH. 


L'analyse des deux approches fournit le Tableau 2 de synthèse suivant: 


Tableau 2 : Comparaison des résultats des deux approches 


STRUCTURE FOUNDATIONS er 


Section | Section Lengh | Widt | Dept 
Column | Beams t h h (DH) 
(mm) | (mm | (mm 
) 


ie nel 464.64 | 181.96 | 742.5 |352.4 |2000.0 50 
Approach 

Genetic 

algorithm 146.30 | 596.97 | 283.33 | 2000.0 40 
Approach 0 


L'apport de la méthodologie des algorithmes génétiques peut être représentée parle 





rapport des coûts engendrés par chacune des deux méthodes : 
50—40 
R=—=p ~ = 0.196 soit environ un gain de 20% 


Il apparaît donc que cette nouvelle approche est très prometteuse et doit être encore 
affinée par un certain nombre de mesures avant la validation des résultats. 


En effet, pour des raisons d’objectivité, l’ordre de grandeur du gain de coût 
pour cet exemple ne peut étre généralisé pour la simple raison que notre outilest en sa 
premiere version et il n'intègre donc pas toutes les contraintes d'optimisation de 
l’Eurocode en termes de calcul de résistance, de stabilité et de dimensionnement des 
constructions. 


En version finale, un vaste programme de simulations numériques sur différentes 
conceptions peut permettre d'élaborer des « règles de conception » qui aident, en phase 


d'avant-projet, les concepteurs vers des solutions globales optimisées. 


6.Application des AG à une poutre-console en béton armé à inertie 
constante 





Afin de mieux évaluer la pertinence de notre approche numérique, nous avons cherché 
à valider notre modèle sur des cas de référence. C’est pourquoi, cette étude porte sur la 
minimisation de la flèche d’une poutre en porte-à-faux soumise a son propre poids. La 
Figure 5 suivant d’écrit notre problème. Dans ce problème académique, n’ayant pour 
unique vocation que de servir de validation pour les développements de codes numériques 
d’ optimisation, nous cherchons à optimiser la forme de la poutre pour minimiser la flèche. 
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Figure 5 : Schéma de la poutre en porte-à-faux 
Cette étude est découpée en deux parties. Dans un premier temps, nous considérons que 


la section de la poutre peut prendre uniquement deux valeurs distinctes et nous 
chercherons la position optimale du changement de section qui minimise la flèche. Nous 
proposons un développement analytique avec une validation expérimentale qui nous 
permettra d'optimiser la forme de la poutre. Nous validerons sur ce cas notre modèle 
numérique, basé sur un algorithme génétique. Puis, nous utilisons notre code pour trouver 
la forme optimale de la poutre dont la section peut maintenant prendre des valeurs dans 


une plage donnée. Dans ce dernier cas, la résolution analytique n’est plus possible. Ces 
deux études peuvent servir de cas de référence pour valider l’approche numérique pour 
l’optimisation automatique de structure. 


Dans ce cas d’étude, l’équation de la déformée de la poutre de longueur Í encastrée à 
une de ses extrémités et libre à l’autre, peut être classiquement décrite par les équations 
de la résistance de matériaux. Le changement de tronçon est identifié par le changement 
de section de la poutre, entre Siet S2, positionné à la distance a de l’encastrement comme 
présenté sur la figure. Les équations de la déformées fiet f2 des deux tronçons de la poutre 
[O,a] et [a,l] sont déduites de la double intégration des équations des moments fléchissant 
de ces deux même tronçons. L’objectif de l’optimisation de cette poutre est de trouver la 
valeur de à qui permet d’observer le déplacement le plus faible de l'extrémité libre de la 
poutre. Ainsi, la validation des approches numériques sur ce cas, portera sur l’estimation 
de ce paramètre a. 


Min { f2(x,a) : a [0,L] |x= L} 


Sur la Figure 6, nous avons tracé la courbe théorique, issue de cette approche analytique, 
ainsi que les valeurs expérimentales de la flèche en fonction de la valeur de a. Ces valeurs 
expérimentales proviennent de mesures réalisées sur des clichés photographiques de la 
déformée d'une poutre en plexiglas. On observe dans ces deux cas que la flèche admet 
bien un minimum pour une valeur particulière de a. Nous évaluons cette position du 
passage entre la grande et petite section de lapoutre à une distance de a = 0,446*/ de 
l'encastrement, ceci dans le cadre d'une poutre en plexiglas de masse volumique 
p=1250kg/m3 et de module d'Young E= 2,2 GPa. Nous avons confronté notre approche 
numérique à ce cas d'étude. L'estimationdu paramètre a optimal ainsi obtenue est égale à 
0,45*/7, soit un écart inférieur à 1% par rapport à la valeur de référence. Nous pouvons 
donc considérer que notremodèle numérique est validé sur ce cas académique. 
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Figure 6 : Evolution de la flèche en fonction du paramètre a 


Notre approche numérique étant validée sur le cas d'étude 1, nous avons souhaité obtenir 
la forme optimale de la poutre dans le cas, cette fois ci, ou la section peut varier dans 
l'intervalle fixé [S1; S2]. Nous observons que la poutre optimale dans ce cas voit une 
variation de section régulière entre S1 et S2 située autour de la valeur précédente de a= 
0,45*/, voir la Figure 7. La section de la poutre est ainsi égale àS1 lorsque l'on se situe a 
une distance inférieure a 0,3 m de l'encastrement et égalea S2 quand nous regardons la 
poutre a une distance supérieure a 0,6 m. 
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Figure 7 : Variation de la section S1 et S2 
Cette étude nous a permis de trouver la forme optimale minimisant la flèche d'une poutre 


soumise à son propre poids. Nous avons étudié deux cas, définis par des contraintes 
différentes sur les valeurs des sections possibles, qui permettent d'obtenir une référence 
supplémentaire en vue de la validation des outils numériques d'optimisation. 


Comme dans tous les pays, les ingénieurs génie civil au Maroc font recours a des méthodes 
d'optimisation pour dimensionner les structures métalliques. A ce titre, nous proposons une 
méthode d’optimisation du poids des structures métalliques formées de poutres en respectant les 
contraintes de dimensionnement définies par la réglementation Eurocode 3. Cette méthode 
d'optimisation, basée sur l'application des Algorithmes Génétiques, consiste à minimiser une 
fonction objective représentant le poids de la structure et dont les variables de conception sont les 
sections des poutres. Ces sections sont choisies à partir d’un catalogue de profilés. Les contraintes 
de dimensionnement englobent la résistance mécanique des sections et des éléments, la limitation 
des déplacements aux nœuds et la limitation de la première fréquence propre de la structure. 

> Dimensionnement optimal d’un portique plan 


Dans cet exemple, nous présentons, les résultats d’optimisation par les AG du poids d’un portique plan de la figure 
1. Les sections des poteaux sont de type HEB et celles des traverses sont de type IPE. Les poutres qui devront avoir la même 
section font partie du même groupe qui sera par la suite une variable de conception. Ainsi, on aura quatre variables de 
conception (tableau 1). Le matériau utilisé est l’acier S235. Les contraintes d’optimisation sont données par les équations 2 à 
6. Les déplacements maximaux sont Uxmax=1.8cm et Uymax=2.4cm. La première fréquence propre dans la structure doit être 
supérieure a 6Hz. A ces contraintes nous ajoutons des contraintes technologiques dues aux assemblages (fig. 3). Elles sont 
exprimées par: b; < b; 


où b; et bj sont les largeurs des barres 1 et J. 
4 5 
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Tableau 1 : Variables de conception 
Groupe Eléments Type de section 








1 HEB 
§ 2 HEB 
3 IPE 
4 IPE 
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Fig. 1: Portique à 2 étages Fig. 2: Maillage du portique Fig. 3 : Contraintes technologiques 


Tableau 2 : Solution d’optimisation du portique plan Tableau 3 : Variables design 





Varablesde | Variable de conception Eléments 
cect on Eléments Solution optimale 1 123527 
1 1 8 HE200B 2 2-4-6-8 
2 07 HE140B 3 9 à 16 
3 3 4 IPE160 A 
4 5 6 IPE360 Tableau 4 : Paramètres de conception 
lère fréquence propre 6.12 Hz Grandeur valeur 
Poids 1300.56 kg Flèche horizontale maximale 1.2 cm 
Taille de la population : 80 - Génération maximale : 60 Flèche verticale maximale 2.0 cm 
Probabilités de croisement pe =0.80 Fréquence propre minimale 4 Hz 


Probabilités de mutation pm = 0.002 


> Dimensionnement optimal d’un portique tridimensionnel 
Dans ce deuxième exemple nous traitons un portique tridimensionnel à deux étages (fig. 4). Le 
portique est sujet à son poids propres et à des charges nodales (F = 6.25KN) et réparties (qı = 
1.5kN/m et q2 = 6kKN/m). Le matériau utilisé est l’acier S235. Les variables et les paramètres de 
conceptions sont donnés aux tableaux 3 et 4. 























Fig. 4: Portique spatial 
Fig. 5: Maillage du portique 


La solution optimale est donnée au tableau 5. Cette solution correspond à un poids minimal de 
la structure de 3456.28 kg et à une première fréquence propre dans lastructure de 4.035Hz. 


Tableau 5 : Résultats de l’optimisation 


Variable de conception Solution optimale 
1 HE280B 
2 HE160B 
3 IPE220 
4 IPE180 


Uiar +1.35 cm 
Uma -0.98 cm 


Masse = 3456.28 kg 
Fréquence = 4.035 Hz 





A travers ce papier nous avons essayé de montrer que l’utilisation des algorithmes 
génétiques, comme outil d’aide a la conception pour guider l'ingénieur vers des solutions 
optimales, est une approche prometteuse. En effet, ils n’exigent aucune condition sur la fonction 
objective (linéarité, dérivabilité, etc.) et peuvent prendre en compte certains contraintes, portant 
par exemple sur la faisabilité et l’esthétique. Ils sont peu coûteux en temps machine et 
permettent de trouver des solutions potentielles proches de l'optimum. Ce qui justifie leur 
utilisation a la place des méthodes déterministes de type gradient ou de programmation linéaire. 
Cependant, les nombreux paramètres qui les contrôlent sont délicats à régler notamment les 
probabilités de croisement et de mutation peuvent faire varier radicalement la vitesse de 
convergence est la qualité des résultats. 

Il existe de nombreuses méthodes qui permettent de résoudre un problème d’optimisation. 
Ces méthodes peuvent se scinder en deux groupes : 

- les méthodes dites déterministes (méthodes exactes). 
- les méthodes dites non-déterministes (méthodes approchées ou méthodes stochastiques). 

Les méthodes déterministes, comme celle du gradient, sont propices à une recherche 
d’optimum local mais ne permettent pas de sortir des creux (puits) pour rechercher un optimum 
global. 

Les méthodes non-déterministes, comme celle de Monte-Carlo, tentent d’éviter une 
convergence vers un optimum local. Cependant, elles ne permettent pas une recherche locale 
efficace. De surcroît, la plupart de ces méthodes ne permettent pas facilement une optimisation 
de problèmes dépendant de variables discrètes. 

De ce fait, les algorithmes génétiques (AG) s'avèrent les mieux adaptés puisqu'ils offrent de 
nombreux avantages par rapport aux méthodes classiques d’optimisation. 

Cette étude nous a permis de proposer une formulation du problème d'optimisation globale des 
structures du Génie civil en général et la construction métallique en particulier. 


Perspectives 


Certes, la méthode présentée dans ce rapport a donné des résultats remarquables et pertinents mais la recherche 
continue ! 

En ce qui concerne les perspectives, il nous semble qu’il reste encore de nombreuses possibilités 
d’exploration, notamment, pour permettre la mise en place de critéres de convergence réellement utilisables. Cette 
notion de critère de convergence nous semble primordiale, c’est sans doute l’aspect le plus limitantactuellement 
pour la diffusion des AG. 


Nous prévoyons élargir l’application de cette nouvelle méthode aux structuresquelle que soit la nature des 
matériaux utilisés (béton ou métal) et dela développer en utilisant des interfaces graphiques mieux améliorées afin 
de faciliter la tâche aux utilisateurs. 


Nous pouvons envisager deux principales directions, l’une portant sur l’étude des propriétés des paramètres et 
leurs rapports aux algorithmes d’optimisation, et l’autre sur la conception de nouvelles techniques d’ajustement 


automatique de ces paramètres. 
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